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Unimodular gravity with external sources. 
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Abstract: The only allowed source of the gravitational field in the unimodular 
£j theory, invariant under area-preserving (transverse) diffeomorphisms as well as Weyl 

^ transformations, is just the traceless piece of the energy- momentum tensor. This 

fact notwithstanding, the free energy produced by arbitrary sources (not only static 
ones) is identical to the one predicted by general relativity. This encompasses all 
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weak field tests of gravitation. 
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1. Introduction 



The purpose of the present paper is to study the basic unimodular gravitational 
action (U from now on) which in the Einstein frame corresponds to a unimodular 
metric, with determinant 

9e = det 9% = 1 (1.1) 

(in Lorentzian signature we always write g to mean the positive number \g\). The 
coupling to matter fields has been worked out in the second reference of [1] using 
Weyl transverse symmetry (WTdiff) as a guiding principle 2 . Its main characteristic 
is that the potential energy does not couple directly to the gravitational field. In 
general relativity (GR) this coupling is through the square root of the determinant 
of the metric, which in our case is set to unity in Einstein's frame. 

This seems promising as a solution to the direct gravitational constant problem, 
namely why vacuum energy seems to weigh much less that other forms of energy as 
if the equivalence principle did not hold for it. This is logically independent of the 
question as to why the universe is expanding in an accelerated way. Nevertheless, as 
we shall see in detail in the following, the interaction of arbitrary external sources is 
exactly the same as in GR. 

The simplest action, using a scalar field as an example, reads 



S = j d n x (-M n ~ 2 R E + \gfd^d^ - Vttfj 



d n x{ g™ 



v{4>) 



= S g + S m (1.2) 
In the second formula we have written the action in an arbitrary frame, using 

9%, = 9i» (1-3) 

It is to be stressed that this is not a field redefinition, because it is not invertible. 
Actually, 

&S& = MtfSgafi = g~» Q + 5X) - Vv) (1.4) 

It is not possible to recover the generic metric out of Einstein's unimodular metric. 
The action is then invariant under local Weyl transformations 

~9,u = e 2 °^g, w (1.5) 



2 This theory was first discovered in a flat space setting in [2] in the context of a general anal- 
ysis of TDiff invariant theories. Those are invariant under transverse (area preserving) diffcomor- 
phisms only. The first such theory is due to Einstein himself [10]. Other relevant references are 
[6][15][16][11][13]. 
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as well as under area preserving diffeomorphisms (id est, those with unit jacobian 
determinant). This is precisely the symmetry that was dubbed WTDiff in previous 
papers. 

The gravitational part of the equations of motion (EM) can be written in the 
manifestly traceless form (which is much more convenient than the one in [1]). In- 
stead of deriving them from the variational principle (which we also did) it is much 
better to use the fact that the Euler-Lagrange equations behave as forms in con- 
figuration space [4]. The Jacobian factor just instructs us to take the trace after 
performing the point transformation on Einstein's equations. This yields in general 

R*p ~ ~Rg a p + (n - 2) (V a V p a - V a aV p a) + (VV - (Va) 2 ) g a p = (1.6) 

n n 

Using now the fact that in our case 

e 2a = g^ (1.7) 
and including also the contribution of the matter field the final result is gotten 



1 (n-2)(2n-l) ( V»gV v g 1 (Vff) 2 \ n-2( V ti V v g lV 2 g 

n An z \ g n 9 J 2n \ g n g 

M 2 ~ n ( 1 \ 

+— W - - AW 9m) = M '~ n { J % + -C) 

It can ve checked that they are Weyl invariant as well. We have included in the source, 
besides the matter part, J™, a gravitational piece J* coming from the space-time 
dependence of the metric determinant. This means that we are defining the sources 
in the matter part 

9* -C = (i-9) 

as well as in the gravitational piece 

g~ (r»» -\R9r + M2 ~ n J9 m) = ( L1 °) 

Taking the contracted covariant derivative and using the contracted Bianchi identity 
corresponding to the Weyl-transformed Einstein tensor 

= v M far - \v v ) = o (i.ii) 

yields the gravitational piece of the Ward identity. 

V„ ( R^ - -Rg" v \ = — -V P R (1.12) 
V n J 2n 
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In this work we shall stick to the standard definition of energy-momentum tensor 
in general relativity (GR) which for an scalar field with minimal coupling means 

T, v = d^d v <t> - {^-g^d^d^ - V(<f>)^j 9ia> (1.13) 

Its trace is 

T ee f^Taf, = nV- T ^ r (V0) 2 (1.14) 
so that the piece of the source that depends on the scalar field is precisely 

JZ = T ilv --T 9llv = d^d v <t> - - g ap d a 4>dp4> g, v (1.15) 



n n 



and does not include the potential energy. It seems plain that 

y frpfMU _ 1 T = l_ V v T , x 16 n 

\ n J n 

We shall see in a moment that this is so in general precisely owing to the Ward 
identities of the area-preserving diffeomorphisms plus Weyl symmetry. To be specific, 
Ward identities guarantee that there is a function 6 such that 

We have to be careful however. The variation of the action under a variation of 
the scalar field reads 

5S m = J d n xg^g^dM - V'b<\> = - J d n x (d v (g* g^d^) + V') <50 = 
It follows that the EM of the scalar field have changed; now they read 

W + ,-iV = ^^» (i.i9) 

2n g 

It is then not difficult to check that 

V, [g^ (v>V>-^(V0)V^) =^V (^+^(W)j=ff^VT 

(1.20) 

which modulo EM implies the desired result. Indeed, for any arbitary constant value 

^ x B= {v ^t + ^t (L21) 

the extra term is multiplied by the EM. 
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In this paper we would like to restrict our attention to the quadratic action 
corrresponding to fluctuations around a given flat background coupled to arbitrary 
external sources. A detailed analysis of the equations of motion will be performed 
in order to compare them to those of GR. The main conclusion of our work is that, 
although the only allowed source of the gravitational field in the unimodular theory 
is just the traceless piece of the energy-momentum tensor 

J av = Tau Tg^v (1.22) 

n 

the full free energy produced by arbitrary sources (not only static ones) is identical 
to the one predicted by general relativity. This encompasses all weak field tests of 
gravitation. 

Then we analyze the corresponding equations of motion corresponding to pertur- 
bations around an arbitrary background. We are able to show that the unimodular 
ones are a subset of the ones corresponding to general relativity, with the source re- 
stricted to its traceless piece and with vanishing cosmological constant. This refines 
a previous analysis by one of us [1]. There it was argued that Bianchi identities on 
the unimodular theory just imply a first integral of the equations of motion which, 
once used, made the unimodular theory fully equivalent to GR with an arbitrary 
cosmological constant. Here we claim a stronger result, namely that at the level of 
external sources, this cosmological constant must vanish. 



2. Ward identities. 

Let us first review the steps necessary to derive Bianchi identities from the field 
theory viewpoint [1]. We perform a transformation of symmetry in the action 

SS= f d n x^-g5g a p (2.1) 

This should vanish identically whenever the variation of the metric is generated by 
a diffeomorphism 

Sg a fi = V a £/3 + V^£a (2.2) 

that is 

= / d?Xy/g (Vofc + V^„) = -2 / d n x^ a V a jL/ g (2.3) 

J vMI 9 a P J \J\g\°9aP 

The fact that the vector is arbitrary conveys the fact that 

V a /^ = (2.4) 

any covariant lagrangian leads to such an identity, which in the case of Einstein- 
Hilbert is the contracted Bianchi identity 
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V a (V 3 - l -R g a ^j = (2.5) 
Let us now impose transversality on the generator 

V a C = (2.6) 

It is useful to define a one form 

6 = Z^dx" (2.7) 

as well as a two form f2 2 

^2 = l^^u dx" A dx v (2.8) 

The transversality condition now reads 

6 = -25fi 2 (2.9) 

where the codifferential [?] is the adjoint operator of the exterior derivative. Acting 
on two- forms 

5 = *^ 1 rf* (2.10) 

so that its components obey 

(SCh) p = ~^ v (^) vp (2.11) 



The number of independent components of a two form is clearly 

(2.12) 



n 
2 

But the codifferential is nilpotent 5 2 = 0, so that we have to withdraw the three 
forms 

n 2 = 5n 3 (2.13) 

and from them we have got to withdraw the four forms, etc. The final counting of 
undependent gauge parameters is: 



2-3-1 



(...) =n-l (2.14) 



where the relationship J2j (~^) J (™) = has been used. 
Then ( taking into account that for antisymmetric tensors Vt al3 = —fl^ a 
V a V^ = 0) 

5S 



= [d n X (V a & + V^a) = I d n X (V a V'fy„ + V/jV'ftap, - 



'a/3 



/-i r c 
d n xJ\g~\ n pfj V p V a —== (2.15) 
\/\g\ og a p 
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Assuming that Q plB is arbitrary (and as we have just seen, they are either arbitrary 
or else vanishing) we deduce that 

V9 og a p yf\g\ bg ap 

that is for the vector 

e0 = V a ±^- (2.17) 



the condition 



which can be integrated as 



where 



V pq/3 = v /3 p (2.18) 



9 p = V p $ + 7p (2.19) 



7 = lp dx p (2.20) 



is an harmonic form. The number of independent such things depends on the topol- 
ogy of the manifold, and is referred to as the first betti number, bi(M), the dimension 
of the first cohomology group, i7 1 (M). 

In case there are not harmonic forms in the space-time manifold, (which happens, 
in particular, if it is diffeomorphic to W n ), this shows that the Bianchi identity is 
modified in the sense that 

V p@/3 = v /3@p = V ^V P $ (2.21) 

id est, still hold when integrated over the whole of space-time with the Diff 
invariant measure. The Weyl invariance of the action means that 

= I d n x w(x) g a ? (2.22) 

which conveys the fact that barring topological subtleties the trace of the EM must 
be a total derivative 

g«e ^- = d p W (2.23) 
og a p 

In order to get Ward identities out of symmetries of the action, we also need the 
parameters to be independent. The standard method starts with a change of variables 
in the path integral expressing the expectation value of a certain monomial of fields, 
X [g pu ,ifji\. where ipi is a generic representation of matter fields. 

Z (0+|X [g, A] |0_> = e iW = J Vg pv Vij; X[g^ ^] e iS ^ +iS ^^ (2.24) 
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namely 

9pu -»■ Qnv + V M £„ + 
This leads easily to the Ward identity 



(2.25) 



SX [g, ipi 



5S 



(2.26) 

In the particular case X = 1, it states that the expectation value of the classical 
identity should vanish. There may be quantum corrections to the naive identities, 
either in the form of anomalies [14] or even limit cycles [12]- 

3. Free energy with external sources in a flat background. 

Let us start by consider perturbations around a flat background 

9fxv = V»v + K-h^ (3.1) 

up to quadratic order. In momentum space the kinetic term of the unimodular theory 
[2] reads 



SS 
5g ap (x) 



K u _ * u2 



pupa 



o k {VppVva + VpaVup) Q {kuka^pp + f]pak v kp + k u kpT]^ cr + k^k^Tj^p) + 



1 . , n + 2 2 

ylp.vkpk a + Tjp a kpk u ) VpuVpa 



(3.2) 



We shall compare all the time with the GR template which at this order corresponds 
to the Fierz-Pauli (FP) spin two theory. Its kinetic energy piece reads 

8K pi/pa = k 2 {j] pp n va + rj^rj"" - 2if v rj pff ) 

- {Wrf 7 + k v k a jf p + k p k a v up + k u k p r] pa - 2k p k u v pa - 2k p k c rf v ) 

This can be easily expressed in terms of the Barnes-Rivers projectors (reviewed 
in the appendix) 

k 2 , „ _ k 2 



K u 

e 

2n 



(2P 2 + 2P S + 2Pi + 2Pq 



(2Pi + 4P " ) + 



{V^TPo* + 2P -) 



n + 2 
An 2 



k 2 {-P 2 - 



n-2 n 2 -3n + 2 
4^2 ^ ' 4n 2 ^° 



k 2 ((n - 1) P S + Vn~^lP x + P?) 
n-2 



An 2 



'V^lPn 



(3.3) 



(The well-known ADM [5] formalism could be used just as well). It is plain that 

(3.4) 



pupa 

£.k = => K^ p k a = 
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Following [2] , we are going to gauge fix WTDiff by adding a gauge fixing 

L gf = V K^ pa h pa (3.5) 

where 

K 9f = ^Pi + ^ ((n - 1)P S + P«T + ^1P X ) (3.6) 
This corresponds in position space to a gauge fixing 

^/ = ^jPa = - °^) 2 ( 3 - 7 ) 

(where A is an arbitrary mass scale) which is admissible, because it can be reached 
uniquely through area-preserving diffeomorphisms 

d a F a = 

5F a = -n 2 ^ (3.8) 

The ghost system gets complicated owing to the fact that the gauge parameters 
are not all independent. For a full analysis confer [9]. It is however irrelevant for the 
purposes at hand, which are purely tree level. 

At the end of the day we are left up with 

4< T = k'P 2 + ((„ - l)m' - n -^e} PS + (m 2 - " 2 ",^" + ^ 2 ) K + 

The euclidean propagator is then given in this gauge by 
, 9a M 4 1 f / , n 2 -3n + 2 j9 \ , /. , „ n-2,.. 



m 



' U *k 2 )V^lP* \ (3.10) 



n 2 



Any coupling of the gravitational fluctuation to an external source 

S in t= I d n xJ^ u (3.11) 



has to comply with 

= 5S int = J d n x J„ v (d^QP + d p d p % + u{x) V) (3.12) 
which is only possible when the source obeys 

= o 

9 M J^ = (3.13) 
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This means that it should be related to some conserved symmetric tensor (which a 
priori could be different from the usual conserved energy-momentum tensor although 
we shall prove it to be the same) by 

1 



J, 



in/ 



T, 



I I III!' 

n 



(3.14) 



The free energy (or effective action) after the gaussian functional integration just 
reads 

1 



W[J] = \f d n xd n yj; u {x) A' 



(3.15) 



The easily proved identities 



tr JP 1 J = 



2) 



J pa = T, 



1 



hi/ 



tr (JP 2 J) = |T^| 2 - 



n — 1 
1 



n — 1 



OuuT 
T 1 2 



convey the fact that 
1 



W[J] = -(27T 



2/i 



where 



C(k) 



1 



This yields the free energy 

1 



W[T] 
1 



2 

(2tt) 2 " 



(27T) 



2/i 



d H k 



(n - l)(n - 2)k 2 
1 /, , „ 2 



(3.16) 

(3.17) 
(3.18) 



A- 2 



J, 



/UW I 



d n k 



k 2 



T, 



[IV I 



n 



n(n — 2) 



|T(fc)| : 



\T(k)\'- 



(3.19) 



This result is already implicitly contained in [2]. This free energy (when expressed 
in the second form) is exactly the same as the prediction of General Relativity. It 
is only in the nonlinear regime that some differences netween both theories could 
possibly be found. 



4. Free energy with external sources in a general background 

In the presence of a vacuum energy flat space is not an allowed classical solution to 
perturb about. Let us now consider instead an arbitrary background g pu . The action 
reads 

(n-l)(n-2) (Vs) 2 ' 



S { 



u 



d n x g™R 



An 2 



g 2 



(4.1) 
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The lagrangian reads 



- i D . (n-l)(n-2)(V</) : 



4n 2 



4 2 



i? + k -i? - /i Q/3 i? a/3 - V 2 h + V a V«/i a/3 + 



n 



2ht>" (-V 2 h vfi + v^v, v + v p v^ v - V^/i) 

-2V a hV a h + 2V a hVph aP - 2V a h ai/ V l/ h arT + V a h av V a h au + Ah p a h va R uP 

4 
— / 

n 

(n-l)(n-2) f(Vg) 2 



-ft (/^i^ + V 2 /i - V Q V^ a/3 ) + 

I 



4n 2 



r 



+ re 2V u /i 



2 ( 


^ 2 


n V 


n 












h al3 h a p | i? 



+ 



9 9 



+ 



K 



! ^-2/i^V^^ + (V/i) 2 - 2(W M /i + h^V^K^ ^ " g g 

The only way the linear term can vanish 3 is by restricting either the allowed 



(4.2) 



fluctuations or else the allowed backgrounds through 



-R - h^R aP - V 2 h + V„V«/^ + 2V„/^ - h^y^ Vu ~ 9 
n 



V 2 h + V a Vph al3 + 2V„h- 



= 



9 



9 9 



For maximally symmetric backgrounds 



R 



2X 



n-2 



9txv 



the condition reads 



V 2 h - V a Vph aP + 2V ll h 



= o 



(4.6) 



(4.7) 



9 9 9 

A simple solution consists in restricting the background to be unimodular by itself, 
id est, 

9 = 1 (4-9) 



3 Let us dispose of a fine point. Usually the integral of a covariant derivative vanishes because it 
can be written as 



1 



(vW) 



so that 



J d n x^V„V = J <Txdp (y/gV") = 



(4.3) 
(4.4) 



assuming vanishing physical effects at the boundary. This is not true anymore with the unimodular 
measure. What can be written instead is 



J d n xgi 



n-2 



(4.5) 
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in which case the offending terms again either vanish or else behave as total deriva- 
tives. 

To summarize, once a unimodular background is chosen the linear term is just 
the equation of motion for the background field 

h^^-^Rg^ =0 (4.10) 

which is the vacuum field equation of the total non-linear background unimodular 
theory. The Bianchi identities force the scalar curvature to be constant. 

We are finally ready for our analysis of the EM of both theories around arbitrary 
backgrounds g^ v and g^ u . Remember that the lagrangians for the U theory and for 
the linearized GR with cosmological constant 4 



L v = V^l^hVph - -V„W"^ + Jv.FV^;- (4.11) 



An 2 n 

Lgrx = \^hV,h - ^V„W'/# + l -V^ p V v h v - (4.12) 
- \v,h^Wh up - R vP hih va + l -hk a ^ - - hi^K^j 

Assuming both backgrounds to be of the maximally symmetric spaces R^ u = 
■^z^g^u, 

Lu = ^^V^hV^h - iv,W"^ + ^FV^; - \v p h» p V»h vp (4.13) 



Lgrx = \^hV,h - \v ,hWh» + l -V»W»V v h v p - Ul p h^Wh up - ^ (j - h af) h a - 

(4.14) 

Sources for both theories can be introduced in the usual way by a linear coupling. 
In the case of GR, the source is just the usual symmetric energy- momentum tensor 
T^ u while for the unimodular theory is its traceless source J^ u which obeys V M J pv = 
V V T. 

The equations of motion of the unimodular theory, dubbed EMU, read 



4 In the full nonlinear theory the cosmological constant is included in an arbitrary energy mo- 
mentum tensor. In the linear approximation this is not the case. 



n 4- 2 - 1- 1-- 
EMU = ~^9^ 2 h - -V 2 V - -VaVph^g^ - 

l^^uh + \v^ a K + ^^aK = J ^ ( 4 - 15 ) 

whereas the equations of motion of general relativity, EMGR are 

EMGR = i V 2 %„ - ^V 2 V - l^a^ph^g^ - \ V,V v h + 

+ ^V„V a /i° = A - v) + ^ (4-16) 

At this point we should remember the result advertised in the introduction, on 
the equivalence of the unimodular theory with GR with an undetermined cosmologi- 
cal constant. We have already analyzed in the previous section fluctuations around a 
flat background and found full equivalence with GR with vanishing cosmological con- 
stant. In order to make sure that ths result is not an artifact of the flat background, 
it is worth to repeat the analysis in this more general setting. 

Two first integrals can be easily obtained. The first one, Iqr, by taking the trace 
of the EMGR: 



I GR ee V 2 h - V a Vph^ -\h= — —^T (4.17) 

n — 2 

(4.18) 

whereas the second one stems from taking the covariant divergence of the EMU 



Uj 



V„ V a V B h a ? --V 2 h+-T 



2n \ n n 



V a Vsh alJ - -V 2 h +-T = T (4.19) 

n n 

where Y is an arbitrary constant. Now assume that the background metric is the 
same for both theories. Since we are trying to check if they are equivalent, this is a 
reasonable ansatz, so we set g^ u = g^ u . After that, we look for a field redefinition of 
the form h^ u = H^ u + aHg^ u that would take one theory into the other , with the 
possible addition of terms proportional to the first integrals which are zero by the 
use of the equations of the motion. This is equivalent to search for constants a, C\ 
C 2 and T such that 

EMGR (H^ + aHg^) + C 2 I GR (H^ + aHg^) = EMU (H^) + CJu (H^) (4.20) 
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V 2 Hg, v - ^ - l - - a (| - l) - C 2 (l + na - a)) + V a VpH a % v (c ± + C 2 + \ - ^ + 

V^V.tf Q - ^ - a (l - |) ) + ff^A Q + a Q - l) + C 2 (l + na)j + 

The system of equations obtained by demanding every factor to be zero is only 
compatible if the cosmological constant A vanishes. In that case, the solution of the 
system is simply 



1 

a = 

n 

2n 



To summarize, 



EMGR [h, v - l -Hg, u y^I GR (h, u - l -Hg^ = EMU(H, U ) |^ =T _i r ^ 

(4.23) 

The physical meaning of what we have proved is that the unimodular EMU are a 
consequence of EMGR when A = only; actually it is the subsector corresponding 
to 

tf? = C - h u g, u (4.24) 

We insist that this is not a field redefinition; is a truncation of GR such that h GR = 0. 
There is no way to build the inverse map from EMU to EMGR. Given the fact that 



h = (4.25) 

is a (partial) algebraic gauge fixing (which does not need ghosts), this shows that 
at the level of the EM the unimodular theory is a truncation of GR with vanishing 
cosmological constant, and with the source reduced to the traceless part of the GR 
source. It is perhaps worth remarking that this does not follow necessarily from the 
fact that the lagrangian is so obtained (gauge conditions can only be used after the 
EM are derived). 
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5. Conclusions 



The main conclusion of our work is that, although the only allowed source of the 
gravitational field in the unimodular theory is just the traceless piece of the energy- 
momentum tensor 

J IIV = Tny Tg^y (^-l) 

n 

the EM of the arbitrary quadratic fluctuations in the unimodular theory are equiv- 
alent to the corresponding EM of General Relativity with the full source T^ v . 

This result holds before any gauge fixing. It could be naively thought that this 
result is just a trivial consequence of the fact that 

9%, = 9~™9^ (5.2) 

Actually this a delusion. The EM corresponding to the lagrangian in which a point 
transformation has been preformed are equivalent to the initial ones only [4] if the 
transformation is invertible, which means that the jacobian must be nonvanishing, 
which is not the case. This means that the correct way of looking at the unimodular 
theory is as a presumably consistent truncation of general relativity, inequivalent to 
it, and one that implies that the cosmo logical constant must vanish. In some sense 
this is not so different as the way superstrings can be understood as a GSO projection 
of the NSR string. 

In the particular case that we are interested in fluctuations with respect to a flat 
background we were able to prove a stronger result namely, that the full free energy 
produced by arbitrary sources (not only static ones) is identical to the one predicted 
by general relativity. 

An important question is whether this truncation will survive quantum correc- 
tions. We hope to report of this in the near future. 
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A. Barnes-Rivers projectors in momentum space. 



Let us briefly state our notation (the same as in [2]). We start with the longitudinal 
and transverse projectors 



kryk 



'a/3 = Va/3 



k 2 



They obey 



e + w = ^ + ^ = ^ = i 

u 2 = u%j} =ul = u (A.2) 



as well as 



tr 9 = n — 1 

tr u = l (A.3) 



The four-indices projectors are 



1 



— i " u,v" pa 

n — 1 

-psw — 

Vn-1 



1 



P ™ = ^==0^ (A.4) 



They obey 



pa pb 

i 3 


SS ab P b 


pa pbc _ 

i 3 


= d ,,<)■''■{>]■' 


pab pc _ 

i 3 


-- (Si/" !";- 


pab pcd 

i 3 





(A.5) 
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as well as 



tr P 2 = *T (^W = 
tr P S = e pa 
tr = uj pa 
tr Pi = 

tr P™ = Vr^T u pa 

f r pws _ \ n 

\Jn-\ 

P2 + P1 + p w + Po = \ (« + 5°s;) (a.6) 

Any symmetric operator can be written as 

K = a 2 P 2 + aiPi + a w P™ + a s P s + a x P X (A.7) 
(where P X = P WS + P$ w ). Then 

K- 1 = -P 2 + -P x + ^P<T + ^ 2 P * ^T^o* (A.8) 

Cfc 2 Qj~\_ (Xi gGjyj Qj y/ CL g Qj'fJJ Oj y/ (Xi g G/UJ Oj y/ 

Sometimes the action of those projectors on tracefree tensors is needed. Defining 
the trecefree projector 

(Ptr) pa XS = \ (S x p 5i + 5X) - \ pu ^ (A.9) 



It is a fact that 



Wp, P ° (Ptr) pa XS = P 2 

n 



1 V^i 



ps p ps ps _^ psw) 

r Q r tr — -<o r Q r Q 

n n 



pw p pw _v puis pw 

r Q r tr — r Q „ „ M) 



1 pw, _ I 

n n 

P\Ptr — P\ 

psw p psw 1 p WS ^ ptu 

p^Ptr = pr - ^Hlp - - — p s (a.io) 
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